home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Collection of Tools & Utilities
/
Collection of Tools and Utilities.iso
/
basic
/
chaosexe.zip
/
XHENON.TRU
< prev
next >
Wrap
Text File
|
1980-01-01
|
1KB
|
45 lines
!HENON MAP
LIBRARY "SGLIB.TRC"
CLEAR
PRINT" ***HENON MAP***"
PRINT
PRINT"THIS PROGRAM GENERATES THE (X,Y) PHASE DIAGRAM FOR THE HENON MAP."
PRINT " X(N+1)=1-A*X(N)^2 + Y(N)"
PRINT " Y(N+1)=B*X(N)"
PRINT" TWO PARAMETERS ARE REQUIRED, A AND B. IF B=1 THE MAP IS "
PRINT"CONSERVATIVE. IF B < ABS(1) THEN THE MAP IS DISSIPATIVE."
PRINT"TRY A=1.4 AND B=0.3 INITIALLY."
PRINT
!
DIM L(1),M(1)
INPUT PROMPT"INPUT INITIAL X , Y VALUES:":XINT,YINT
INPUT PROMPT"INPUT A,B VALUES:":A,B
INPUT PROMPT"INPUT NUMBER OF THROWAWAY ITERATIONS:":INITNUM
INPUT PROMPT"INPUT NUMBER OF TOTAL ITERATIONS:":FINNUM
!
CALL SETXSCALE(-1.5,1.5)
CALL SETYSCALE(-.5,.5)
CALL SETTEXT("HENON MAP","X","Y")
CALL RESERVELEGEND
!
DATA 0,0
CALL DATAGRAPH(L,M,1,0,"WHITE")
CALL GOTOCANVAS
LET X=XINT
LET Y=YINT
FOR I=INITNUM+1 TO FINNUM+1
LET XNEW=1-A*(X^2)+Y
LET YNEW=B*X
CALL GRAPHPOINT(XNEW,YNEW,1)
LET X=XNEW
LET Y=YNEW
NEXT I
CALL ADDLEGEND("A="&STR$(A)&" B="&STR$(B),0,1,"WHITE")
CALL DRAWLEGEND
get key variable
clear
print"press <esc> key to finish"
END